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Computation of the Heavy-Light Decay Constant with NRQCD 
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Non-relativistic QCD is applied for a lattice computation of the heavy-light meson decay constant in quenched 
approximation at (3 — 6.0. Clear signals are obtained for the ground state at large times in the correlators, 
allowing a reliable extraction of the decay constant. Estimating the current renormalization factor by the tadpole 
improvement procedure, we find fs = 164(17) MeV for the B meson with a -1 = 2.3 GeV, while an extrapolation 
to the static limit yields /| tatlc = 247(26) MeV. 



1. Introduction 

In recent years large effort has been directed 
toward a lattice QCD calculation of the decay 
constant of heavy-light mesons. Uncertainties in 
the results, however, are still quite large @. In 
the static approximation for heavy quarks, the 
problem manifests in the presence of large noise in 
the correlators, which makes it difficult to extract 
ground state signals. Smearing techniques have 
not resolved the problem, failing to yield results 
independent of the smearing sizeQ. 

In this work we report on a calculation of 
heavy-light decay constants using non-relativistic 
QCD (NRQCD) j|. Our motivation stems from 
the expectation that noise in correlators should 
be reduced for a finite heavy quark massQ 
and hence cleaner results could be obtained in 
NRQCD. In addition NRQCD allows incorpora- 
tion of 1/mg corrections in a systematic way. 
This is an important point since available results 
indicate that 1 /tuq corrections to the static limit 
is fairly large for the B meson Q. 

2. Simulation 

For the heavy quark we employ the standard 
NRQCD action given by§] 
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A 4 -H 



(n) 



Cl 



S-B 

2m,Qa 



Qxt(l) 



where ff« =H = J2j A- 3 A 3 /2m Q a and = 
H — -ff 2 /4, the latter being a modified choice 
in order to stabilize high frequency modes for 



rriQa < 3. The spin-magnetic interaction term 
is included for keeping consistency of the 1/tuq 
expansion. We use the tree-level value c\ = 1 for 
the coefficient. 

Our simulation is carried out with 40 quenched 
configurations on an 16 3 x 48 lattice at /3 = 
6.0. The heavy quark masses used are rriQa = 
1000,10.0,7.0,5.0,4.0 with the n = 1 action and 

5.0. 4.0.3.0.2.5 with n = 2. The value m Q = 
1000 is taken to compare with the results of the 
static approximation. For the light quark we 
used the Wilson action with K = 0.1530, 0.1540 
and 0.1550. The critical hopping parameter is 
K c = 0.15708(2). 

The heavy-light meson decay constant is ex- 
tracted from the correlator of the axial- vector cur- 
rent and the heavy-light meson in the standard 
manner. For the meson operator we use the local 
form and also employ the cube smearing of sizes 
3 3 , 5 3 , 7 3 and 9 3 , fixing gauge configurations to 
the Coulomb gauge to avoid the necessity of in- 
serting gauge link factors. 

3. Results 

3.1. Signal to noise ratio 

In Fig. [l] we show the effective mass of the local- 
local correlation function 



C{t) =< Q(t) l4l5 q(t)q(0) l5 Q(0) > 



(2) 



for mQa=1000 and 5.0. Clearly the signal is 
far better for mQa=5.0 for which we observe a 
plateau beyond t f» 10 — 12. 
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Figure 1 . Effective mass of the local-local correla- 
tion function for rriQa — 1000 (open circles) and 
rags = 5 (filled circles), both with K = 0.1530 
for light quark. 




The improvement of ground state signals for 
smaller values of rriQa can be qualitatively under- 
stood from the estimate of the relative errorBl, 



5C(t) 



where E(Qq) and E(QQ) are binding energies 
of heavy-light and heavy-heavy mesons. For fi- 
nite rtiQa the negative contribution of E(QQ) re- 
duces the value of the exponential slope from that 
in the static limit where E(QQ) vanishes. We 
found that our data for SC(t)/C(t) are consistent 
with this estimate. Typical examples are shown 
in Fig. |^ where solid lines indicate the slope ex- 
pected from the measured values of the binding 
energies and m^. 



Figure 2. Relative error of local- local correla- 
(3) tion function for mQa = 1000 (open circles) and 

mqa = 5 (filled circles), both with K — 0.1530 
for light quark. 
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3.2. Dependence on smearing size 

In Fig. |^ we plot the raw value of the com- 
bination fpyjmp of the decay constant fp and 
the heavy-light meson mass mp extracted from 
fits of the correlators over the interval t m in < t < 
tmin+i for various smearing sizes. For each group 
of data points t m in increases as t m in = 6, 8, 10, 12 
from left to right. We observe that the estimates 
converge to the same value after t ss 10 for all the 
smearing sizes including the case of no smearing. 
This gives us confidence that the asymptotic re- 
gion is really reached at t s» 10 — 12. We also note 
that the magnitude of errors are similar for var- 
ious smearing sizes. In particular reliable results 
can be obtained without smearing in NRQCD. In 
the following analysis we use the cube smearing 
of size 5 3 and extract /p^/mp from a global fit 
over the interval 10 < t < 20. Other choices give 
similar results. 



3.3. Heavy-light decay constant 

In order to obtain fpsjmp as a function of the 
heavy quark mass toqgi we extrapolate the results 
at three values of K for the light quark linearly in 
1 / K to K c for each rriQ . For the axial- vector cur- 
rent renormalization factor Za we take the value 
Za = 0.65 obtained by applying the improve- 
ment procedure with the coupling 5y(l/a)[f| to 
the one- loop result^), disregarding a small ttiq 
dependence of a few % over our range of rriQa. 
The results for fp^mp are plotted in Fig. ||. Cir- 
cles and triangles are for the results obtained with 
the n = 1 and 2 action in (1). We see that they 
yield consistent values for fp^/mp in the region 
of vriQa where both actions can be employed. 
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Figure 4. 1/tiq dependence of fp^/mp after 
extrapolation the limit K = K c for light quark. 



Figure 3. Dependence of fp^/m P on the fitting 
range t m i n < t < t m i n + 4 for various smearing 
sizes at m,Qa = 5.0 and K — 0.1530. 



In order to see the deviation from the heavy 
mass scaling law fpyjmp = const, we fit this data 
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with the form, 

f P V^ = (/pV^?) static (i - — ), (4) 

m Q a 

and obtain 

(./p\Mp) static = 0.163(17), c = 0.60(20). (5) 

With a -1 = 2.3(1) GeV obtained from the p me- 
son mass on the same set of configurations, this 
gives 

/static = 247(26) (^-^ 1 MeV (6) 

for the B meson decay constant in the static ap- 
proximation. 

Substituting rriQa = 1.8 for the b quark mass 
which is estimated from the T mass including the 
heavy quark mass renormalizationQ, we find a 
substantially smaller value for the B meson decay 
constant in NRQCD: 

^^^(^Gev) 37 "^ (7) 

The large coefficient c=0. 60(20) in (^) leads to a 
significant 1/toq correction (~ 30 %) to the static 
approximation, which is consistent with the re- 
sults obtained with propagating heavy quarks [Q. 
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